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Summary 


• The reduced opening angle technique offers a simple way with minimal 
model dependence to measure cosmic ray energies over a broad range with 
out any normalization uncertainties. 

• The emulsion film and CR39 detectors proposed are well established 
techniques and should perform adequately. 

• The analysis method depends on accurate automatic scanning of the CR39 
plates. UAH have developed such a capability. 

• With the proposed geometry energy measurements to ~ 5 TeV/a can be 
made. 

• The expected iron event rate (E>500 GeV/a) is 10 m' 2 day' 1 . 

• The expected energy resolution, from accelerator calibrations at 200 
GeV/a, is -50% to +80%. Since the absolute flux has some sensitivity to the 
assumed power law index it is essential that good energy resolution is 
obtained. 

•The expected charge resolution is ~ 0.3 charge units for the CNO group 
falling to ~ 1 charge unit for the iron group. 

• A suitable event trigger would be a measurable (> 2 pm) deflected heavy 
(Z>2) fragment. 

• One potential background is electromagnetic dissociation that 
predominantly couples to individual protons or alphas. Although the cross- 
sections can be appreciable such events will not pass the event trigger. 


2 


1. INTRODUCTION: 


5 


2. THE REDUCED OPENING ANGLE TECHNIQUE 5 

2.1 The physics of the technique 5 

2.2 The formalism 6 

2.3 Comparison with experimental heavy-ion beam data. 7 

2.4 Experimental data 8 

3. THE PROPOSED DETECTOR 1 0 

4. THE EXPECTED PERFORMANCE 1 1 

4.1 The fragment distribution. 1 1 

4.2 The Charge resolution 1 4 

4.3 The energy resolution and the maximum energies that can be measured. 1 4 

4.4 Event rates. 14 

4.5 Determination of the absolute intensity of the primary cosmic radiation. 15 

4.5.1 The detection efficiency 15 

4.5.2 Fragmentation processes in the overlying atmosphere. 1 6 

4.5.3 The effects of energy resolution and the geomagnetic cutoff momentum. 1 6 

5. CONCLUSIONS 17 

6. APPENDICES 17 

6.1 Appendix A: Electromagnetic dissociation 17 

6.2 Appendix B: Nuclear cross-sections and their energy dependence (E<10.6 GeV/a). 19 

6.3 Appendix C: Some measured fragment distributions 20 


3 



Figure I q 2 distribution of a-fragments for various projectile at several energies 7 

Figure 2 Accelerator calibrations of the opening angle technique - 

Figure 3 q 2 distribution for a fragments from a range of cosmic ray primary particles 

Figure 4 q 2 distribution for a range of incident cosmic rays where the emission angle has ^ 

been scaled according to equation 5 

Figure 5 Schematic diagram of the proposed detector construction 

Figure 1 1 Target factors for a 10.6 GeV/a ,91 Au on three targets 1 


Figure 8 The yield of all fragments u ith Z > 2 as a function of fragment number for a 

,97 Au projectile """ 

Figure 9 The ratio of the protons yield to the number of events as a function of the 

number of released protons for a 197 Au projectile 

Figure 10 Pseudo rapidity distributions for a |;7 Au projectile 


Table 1 Values of the electron scattering radii and associated parameters 

Table 2 Parameters for equation 15 — — ........... I - 

Table 3 Calculated cross-sections for a 500 GeV/a iron cosmic ray incident on a lead 
target. No nuclear structure effects are included. The cross-sections should be taken 

only as a semi-quantitative guide - ^ 

Table 4 The defections in pm at a distance of 0.3 m for a range of species and for y values 

between 500 and 5000 (0.499-4.999 TeV/a) [ 4 

Table 5 The aperture factor for two 1 m 2 detectors as a function of their separation 15 

Table 6 The expected iron event rate for two lm 2 detectors separated by 0.3 m.. 1 5 

Table 7 Measured electromagnetic dissociation cross-sections for 1 Au projectiles and 

deduced values for iron (see text) 

Table 8 Fit parameters for hard sphere models with overlap terms given by equations 2^ 

and 26 7 2C 

Table 9 Iron and lead cross-sections calculated with and without energy dependent 

models 


4 



1. Introduction: 

Accurate measurement of the primary galactic cosmic ray species energy 
dependence in the regime beyond ~ 500 GeV/a is difficult due to the low flux and the 
limitations of energy measurement techniques. However, such observations are essential 
to resolve several questions of current interest such as: Is the enrichment of heavy species 
(Z>6) cosmic rays first reported at higher energies by the proton satellite 1 and then later at 
lower energies 11 real? The results from a previous deployment of the reduced opening 
angle technique are inconclusive' 1 but the authors do point to limitations in the previous 
techniques. Another intriguing puzzle is the energy dependence of silicon cosmic rays. 
Two independent experiments using different experimental techniques indicate that 
silicon is under-abundant. At present the observation is limited by statistics; it could still 
be a three sigma fluctuation. However, if confirmed the current models of acceleration 
and propagation which are species independent are seriously inadequate. 

To progress further the species and energy dependence must be accurately 
measured in a manner that is free from systematic uncertainty. In this report we show that 
the reduced opening angle method offers a simple and relatively inexpensive method to 
answer these questions. 

First we present the physics of the reduced opening angle and indicate the 
expected energy and charge resolution. The proposed detector design is then presented 
followed by the expected performance. Where ever possible simple phenomenological 
expressions that allow ‘back of the envelope’ estimates are given. More details are 
presented in the appendices. The limit of the energy resolution and the expected event 
rates for iron cosmic rays are calculated. Salient points are summarized in the 
conclusions. 

2. The reduced opening angle technique 

2. 1 The physics of the technique 

The method relies on the experimental observation that the fragments from an 
high energy collision, in the rest-frame of the projectile, have a universal isotropic 
Gaussian momentum distribution that is independent of the fragment mass: the decay is 
independent of formation'"’"'. The target nucleus acts only to inject energy into the 
projectile, raising it to a state of high excitation and causing it to ‘explode’ into a number 
of fragments. There is no memory of how the energy was transferred. The relative 
probability for a given fragment is independent of the excitation process. We note in 
passing that this technique shares most features of the CastagnoIi v method 2 used for 
mesons. 


2 The Castagnoli method is derived from the well known relation tanG = y L tan0*/2 
between the lab. angle 0 and the c.m. angle 0*. This is obtained by the approximation that 
6* =sqrt(p* 2 +m 2 J ~ p* 2 in the c.m.s. This approximation is not applicable in the case of 
evaporated fragments from nuclei since p* ~ p F « M N . Hence s* ~ M N 
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2.2 The formalism 

We shall closely follow the formalism of Ichimura et aT. We first restrict the 
discussion to a-particles and then show how this may be extended. 

With the a opening angle technique the primary energy per nucleon (E 0 ) is given 
by: 


E 0 * Po 


Pra JJra) 

sin0 a ~ 0 a 


1 . 


with <P Toc > » 70-90 MeV/c, where 0 K is the emission angle of the a in the laboratory 
system. Though equation (1) is easy to use, the energy being determined from the 
quotient of <P Ta > to 0*. there are several practical problems: 


1) Fluctuation from the njean 

2) The energy is influenced strongly by the most forward a-particle 

3) Not all the fragments are a-particles. 


The third problem is most critical especially for light nuclei: oxygen to silicon. A 
solution to these difficulties has been developed based on the theoretical model of 
Goldhaber". Let P f be the total momentum of a fragment with mass number A* produced 
from the projectile with mass number A,,. Goldhaber showed that the mean square 
momentum <P 2 ,> in the projectile rest frame is given by: 


< P]) = (Pl) 


4.-1 


2 . 

where <P 2 N > is the mean square momentum of a nucleon. It is approximately related'" to 
the Fermi momentum p f by 



180 MeV/C 


3. 


Although equation 2 is derived from simple assumptions it reproduces heavy-ion-beam 
fragment data welT”. Using this equation we immediately obtain the following relation 


(Pi)={Pw) 


04 ,- 1 ) 


4. 


with <P 2 TN > = (2/3) <P 2 n >; the transverse momentum is expressed per nucleon. 

We have so far only considered a-particles. To extend the model to the full range 
of fragment species we introduce a reduced angle 0 given by: 
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6 = 6j 

5. 

The reduced angle calculated from equation 5 is independent of the mass of both the 
projectile and the fragment Hence while the angular distribution 9 depends on the 
projectile, energy as well as the participant masses by using the reduced angle a 
considerable simplification occurs: 

9(E„,e f ; Af Ap)d0 2 f = 9 (q)dq 2 

6 . 

with q = P L y L 0 ~ y L 0, where P L and y L are the velocity (c = 1 ) of the projectile in the lab. 
system and its Lorenz factor. 

This is a crucial result: we can treat equivalently any fragment p, a, Li, Be, ... 
Consequently the fluctuations in the average transverse momentum will be considerable 
reduced from those that come with only measuring individual a-particles. 

2.3 Comparison with experimental heavy-ion beam data. 

The validity of the scaling law (equation 6) can be seen in Fig. 1 . The scaling law 
holds for a range of species over a wide energy range. The shape may be reproduced by 
the summation of two Gaussian curves fit to the experimental data with the constraint that 
t,+t 2 = 1. 


Figure 1 q z distribution of a-fragments for 
various projectile at several energies. 


o’ - stfe* 

<f> = z’,a 2 e ( "° ,I,J) + r 2 tf 2 ,e (-a ' rfI) 

7 . 
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8 . 


with a l 2 - 


M 


N 




r, t, =0.7252, x 2 = 0.2748, a, = 98.5 MeV/c, ct 2 - 239.6 MeV/c. 


u 


We now calculate the mean value of lnq (= lnp L y L 0) 


<q>=\Q\nq<p(q)dq 2 


1 

— / + In 



2 



where y (=0.5772) is the Euler constant, and 
cr 0 = cr*' a 2 2 = 125.8 MeV / c 
so that 


7 


2 -2(ki0) 

I 

0 


= \l 1+ r n e *y Q e 


-(M 


Hence the average value of 9 is given by 


9. 


10. 



Yl 


with 


11 . 


r* =: rr Le ~ r ' 2 = 01422 

m n 

12 . 

The mean value of <ln0> is obtained experimentally by 

(in0>=— 

n / f-i 

13. 

where iif is the number of fragments. 

From a practical point it is important to to exclude wounded protons (protons that 
receive significantly more energy transfer) and or tz^s in the measured data. To reduce 
their number minimum ionizing events with emission angles five times the smallest are 
excluded in addition to a cut made on the charge sum restriction. 


2.4 Experimental data 

This applicability of this technique has been studied using accelerator beams and 
works well as shown in Figure 2. 
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Figure 2 Accelerator calibrations of 
the opening angle technique. 

The standard deviation a is of 
magnitude 0.2 - 0.3 for log 10 E est /E tnie . 
Taking a putative c of 0.25 the 200 
GeV/a data lies between 1 12.3 and 356 
GeV/a. That is a one sigma error of - 
44% to +78%. 


The q 2 distributions for a 
events from a range of incident 
cosmic ray projectiles is shown 
in Figure 3. There is a 
remarkable species and energy 
independence to the distribution. 

Figure 3 q 2 distribution for a 
fragments from a range of 
cosmic ray primary particles. 
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FIG. 7. Resolution of our 
heavy-ion beams at several ei 


eneryy determination calibrated by 
icrfies. 



For the present purposes we are also interested in the q 2 distribution for heavy 
fragments which is shown in Figure 4. The angle scaling works remarkably well clearly 
demonstrating that the technique can be used for all fragment species so long as the 
fragment mass is known. 
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Figure 4 q 2 distribution for a 
range of incident cosmic rays 
where the emission angle has been 
scaled according to equation 5. 


3. The Proposed detector 


\ Incident cosmic ray 



sided 75 micron 
emulsion on 250 
micron support 


The proposed detector is 
relatively simple. It would 
consist of a number of modules 
as shown schematically in 
Figure 5. To reduce 
development and fabrication 
times and to minimize the 
instrument cost each module 
would be the same size (40x50 
cm) as used for JACEE 
experiments. 

UAH has developed 
digital tracking techniques so 
that the cosmic ray species and 
the trajectory can be determined 
automatically by scanning the 
CR39 plates. Using this 
information the expected 
location of the pool of 
fragments on the emulsion is 
examined and the fragment 
charge and angular deviation 
determined. Hence, as discussed 
in 2.2, the energy of the cosmic 
ray can be determined. 
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Figure 5 Schematic diagram of the proposed detector construction 

For the design proposed here lead target foils are suggested. However it would 
also be possible to use silver emulsions as a target. This was the approach adopted by 
Ichimura vl et al. 

4. The expected performance 

4.1 The fragment distribution. 

Nuclear interaction fragment cross-sections have been measured at energies up to 
10.6 GeV/a for ultra heavy nuclei by Garrard et al. vl " and Waddington et al. K . A detailed 
discussion of the total and partial charge changing cross-sections is presented by Nilsen et 
al\ Since the total charge changing cross-section provides an upper bound to the 
measurements and since new values based on nuclear charge radii are presented by Nilsen 
et al. We will first discuss them. 

The conventional" approach to nuclear cross-section calculation has been to take 
the radius" as a constantxA 173 . There are several problems with this approach. Firstly 
hydrogen does not naturally fit into the series and more importantly electron scattering 
form factors indicate that nuclear charge radii do not scale directly as A 173 . As many 
charge radii have been measured a more useful approach is to use these data (Nilsen et 
al x .) In Table 1 we show some measured radii taken from their paper. 



H 

Li 

C 

Al 

Fe 

Cu 

Kr 

r 0 =RVA 173 fin 

1.034 

1.64 

1.36 

1.25 

1.27 

1.28 

1.22 


1.35 

2.57 

3.09 

4.05 

5.16 

5.39 

5.91 

R, fin 

1.034 

3.13 

3.11 

3.75 

4.85 

5.11 

5.34 

R s »=sr, 

1.32 

4.00 

3.97 

4.79 

6.19 

6.52 

6.82 


Ag 

Sn 

Xe 

La 

Ho 

Au 

Pb 

r 0 =R/A 173 fin 

1.23 

1.20 

1.21 

1.21 

1.22 

1.17 

1.19 


6.45 

6.63 

6.86 

6.99 

7.40 

7.85 

7.99 

R,fin 

5.87 

5.90 

6.15 

6.27 

6.69 

6.81 

7.04 

R n = SR e 

7.39 

7.53 

7.85 

8.01 

8.54 

8.70 

8.99 


Table 1 Values of the electron scattering radii and associated parameters. 

Nilsen find the best fit to the energy-independent charge changing cross-sections 3 4 5 6 
is 

a(R p , R t ) = n x [l? p + R t - (320 ± 0.05)] 2 


3 R* is the measured electron scattering radius. 

4 is a conventional radius using a normalization of L35xA 1/3 as used by Westfall. 

5 R* is the nuclear radius computed as a scaling factor (1 .277) times the charge radius. 

6 Cross-sections are usually measured in bams (lO'^m 2 ) whereas nuclear radii are measured in fermi (1(T 15 
m). Hence 1 bam = 100 fin 2 . 
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14 . 

The radii used are the matter radii (RJ taken from Table 1. For iron cosmic rays hitting a 
lead target we expect a total charge changing cross-section of 4509 mb. Nilsen et al. have 
extended their model to incorporate logarithmic energy dependence. However, the 
relatively small energy change from 0.6 GeV/a to 10.6 GeV/a produces a set of very 
different parameters see appendix B. Moreover, Nilsen et al. express reservations (see 
their paper for an explicit discussion) about extrapolating these expressions to higher 
energies. Hence we will use equation 14 to calculate our reference cross-section. 

Estimation of the partial cross-sections is more difficult. Some authors*"’' 1 " 
believe that at high energies in the regime beyond reaction thresholds and resonance’s the 
cross-sections will have only a small logarithmic energy dependence. Here we will 
present some estimates based on Nilsen et al. The models are relatively crude containing 
no nuclear structure information and are unlikely to give a good description of data at the 
high energies which are of primary interest in this report. However, they do provide some 
guidance. Measured fragment yields for 10.6 GeV/a gold projectiles are presented in 
appendix C. These demonstrate the limitations of these simple models. 

We compute the partial charge changing cross-section according to the 
prescription presented by Nilsen et al. Here AZ is the change in the projectile charge. It 
can be positive (a pickup-reaction) although it is predominantly negative (a stripping 
reaction): 


cr*z(A A T ,K,AZ) = p^A 


+ Af - p 2 )(1 + — )|AZT M,+ ^ /ft+ ^ ,/ ' t+ * / ^ 
K 


15 . 


The parameters are given in Table 2. A„ T are the projectile and target masses respectively 
and K is the kinetic energy (Total energy - n^c 2 ) in GeV/a. 


p. 

21.2±0.5 mb 

p 2 

1.08±0.15 

p 3 

(0.485±0.014)A GeV 

p 4 

0.094±0.013 

p 5 

1.11±0.02 

p« 

10.8±1.6 

p 7 

(0.85±0.03)A GeV 


Table 2 Parameters for equation 15. 


Nilsen et al. also presents simple phenomenological algorithms, with a simple 
energy scaling, to calculate the expected fragment distribution yield. No nuclear structure 
information is included and the modeled energy dependence gives asymptotic energy 
independence in the regime of interest here. We also note the cross-section algorithms 
were optimized for the fragmentation of ultra heavy nuclei where one would expect some 
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washing out of the strong structural effects seen in sd shell (Z= 8 to 20). Although the 
cross-sections were derived from relatively low energy (< 10.6 GeV/a) data it is 
anticipated that there will be no serious deviation from the weak energy dependence 
incorporated in the model. However, it should be borne in mind that older attempts to 
extrapolate nuclear cross-sections at lower energies have frequently been found 
inadequate when higher energy data became available. The cross-sections are shown in 
Table 3. 


Cosmic ray energy 

Total cross-section 7 

Total charge changing cross-section 8 

500 GeV/a 

4508 mb 

1869 mb 

Zchange 

Cross-se 

(mb 

ction 

% of the 
charge 
changing 
cross- 
section 9 





% of the 
charge 
changing 
cross-section 

+1 

184 

9.8 

-13 

55 

3.0 

-1 

184 

9.8 

-14 

53 

2.9 

-2 

133 

7.1 

-15 

51 

2.8 

-3 

109 

5.9 

-16 

50 

2.7 

-4 

96 

5.1 

-17 

49 

2.6 

-5 

86 

4.6 

-18 

47 

2.5 

-6 

79 

4.2 

-19 

46 

2.5 

-7 

74 

4.0 

r2o 

45 

2.4 

-8 

69 

3.7 

-21 

44 

2.4 

-9 

66 

3.5 

-22 

43 

2.3 

-10 

62 

3.3 

-23 

42 

2.3 

-11 

60 

3.2 

-24 

42 

2.2 

-12 

57 

3.1 

-25 

41 

2.2 


Table 3 Calculated cross-sections for a 500 GeV/a iron cosmic ray incident on a lead 
target. No nuclear structure effects are included. The cross-sections should be taken 
only as a semi-quantitative guide. 


7 The total cross-section calculated from equation 14. 

* The charge changing cross-section calculated from equation 15. 

9 These are % of the total charge-changing cross-section, (equation 15) which has been greatly extrapolated 
from its know region of validity. NB Since the total hadronic cross-section is calculated from equation 1 4 
to be twice that of equation 15 these values are uncertain to at least a factor of 2. 
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4.2 The Charge resolution 

The charge resolution is determined by the CR 39 etch detectors. Typical values 
of < 0.3 charge units should be achieved for CNO group cosmic rays falling to ~ 1 for 
iron group species. 

4.3 The energy resolution and the maximum energies that can be 
measured. 

The measured energy resolution (one sigma) for projectile energies up to 200 
GeV/a is -44% to +78% as shown in section 2.4. 

For a practical detector the limiting energy resolution is set by the smallest 
deflection that can be measured which is itself controlled by the physical separation 
between the target and emulsion plates (see Figure 5). For comparison purposes we will 
take the following practical conditions: 

1) Target to emulsion separation 0.3 m 

2) Emulsion grain sizes of 0. 1 pm and the minimum deflection that can be measured is 2 
pm. 

Using equation 1 1 we calculate the average emission angle and relate this angle to 
a specific fragment angle using equation 5. In Table 4 we show the calculated deflections 
in pm as a function of energy. If we set our detection criterion as at least one heavy (Z>2) 
fragment we can perhaps hope to measure energies to ~ 5 TeV/a with this geometry. 


m 

H 

Ha 

IB 

mm 

1E8M 

20 Ne 

24 Mg 

mm 

mm 

36 Ar 


500 

85.3 

41.5 

26.3 

22.0 

18.2 

15.4 

13.3 

11.5 

10.0 

8.6 

7.3 

1000 

42.7 

20.7 

13.2 

11.0 

9.10 

7.72 

6.64 

5.75 

4.98 

4.29 

3.6 

1500 

28.4 

13.8 

8.8 

7.3" 

6.1 

5.1 

Esa 

3.8 

3.3 

2.9 

fsm 

2000 

21.3 

10.4 

6.6 

5.5 

QJ|[ 

3.9 

3.3 

2.9 

2.5 

2.1 

1.8 

2500 

17.1 

8.3 

5.3 

Esa 

3.6 

3.1 

fsm 

2.3 

2.0 

1.7 

1.5 

3000 

14.2 

6.9 

Esa 

E 

3.0 

2.6 | 

2.2 

1.9 

1.7 

1.43 

1.2 

3500 

12.2 

5.9 

3.8 

3.2 

2.6 

2.2 

1.9 

1.6 

1.4 

1.2 

1.0 

4000 

10.7 

5.2 

3.3 

2.8 

2.3 

1.9 

1.7 

1.4 

1.3 

1.1 

0.9 

4500 

9.5 

4.6 

2.9 

2.5 

2.0 

1.7 

1.5 

1.3 

1.1 

1.0 

0.8 

5000 

8.5 

4.2 

2.6 

2.2 

1.8 

1.5 

1.3 

1.2 

1.0 

0.9 

0.7 


Table 4 The defections in pm at a distance of 0.3 m for a range of species and for y 
values between 500 and 5000 (0.499-4.999 TeV/a). 

4.4 Event rates. 

The instrument has been modeled as two parallel rectangles placed directly above 
one another and separated by a constant distance. The aperture factor has been 
computed 10 using the techniques of Sullivan”' 1 for arbitrary dimensions. However no 
attempt has been made to allow for the finite detector thickness nor for nuclear interaction 


10 Code is on Smiles in the directory $2$dkal00:[smitha.nuc_int]ap_fac_rectangle.for 
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losses. Hence the numbers are upper limits. A flight instrument would probably use 
JACEE geometry (40x50 cm) but for comparison purposes we here present the aperture 
factors for two square (1 m 2 ) detectors. 


Table 5 The aperture factor for two 1 m 2 detectors as a 
function of their separation. 


To calculate the event rate we will take aperture factor for a sheet separation of 30 
cm, an exposure of 1 day, and calculate 11 the integral iron event rate for a range of 
energies. 


Table 6 The expected iron 
event rate for two lm 2 
detectors separated by 0.3 
m. 


Y 

TeV/a 

Iron event rate/day with 
0.3m between two 
square (1 m 2 ) detectors. 

100 

0.1 

Ill 

500 

0.5 

10 

1000 

1.0 

4 

5000 

5.0 

0.3 

10,000 

10.0 

0.1 

50,000 

50.0 

0.009 

100,000 

100.0 

0.004 


Sheet 

separation 

(cm) 

Aperture 

factor 

(m 2 sr) 

20 

2.17 

30 

1.82 

40 

1.54 

50 

1.30 

100 

0.63 


4.5 Determination of the absolute intensity of the primary cosmic 
radiation. 

The determination of the absolute intensity of the primary cosmic radiation 
depends on a knowledge of three factors (that will be discussed separately): 

1. The detection efficiency 

2. The fragmentation in the overlying air mass 

3. The effects of the energy resolution and the geomagnetic cutoff momentum. 


4.5.1 The detection efficiency 

Often it is difficult to calculate the detection efficiency for a counter experiment. 
One of the advantages of the reduced opening angle technique is the simple geometry. 


11 Code is on Smiles in the directory $2$dkal00:[smitha.bugs.flux]int_flx_rec 
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Hence the predominant uncertainty is the nuclear cross-section. There are several models 
including the one we discuss in appendix B that uses measured nuclear charge radii rather 
than values scaled by the cube root of the average element mass. Ichimura et al ." used 
such a simple scaling model of KaroP v and found that this did give a good reproduction 
of the measured cross sections as determined from the ratio of jets to primary particles. 
Hence the uncertainties introduced by the cross-section calculations are relatively small. 


4.5.2 Fragmentation processes in the overlying atmosphere. 

By necessity any balloon deployment will be below an overlying atmosphere ol 5- 
1 0 g/cm 2 . Since this is the common situation a number of techniques have been developed 
to correct a measured flux to that at the top of the atmosphere. The techniques must also 
account for the flux of secondary particles (from the interaction of primaries with the 
atmosphere) and correct the measured flux accordingly. There are a number of reviews in 
the literature including that employed by Ichimura et al. vl in their deployment ol the 
reduced opening angle technique. Uncertainties from such atmospheric interactions are 
small and quantifiable. 


4.5.3 The effects of energy resolution and the geomagnetic cutoff momentum. 

The factor required to calculate the measured flux from the observed flux depends 
on the geomagnetic cutoff momentum, the energy resolution and the power law index of 
the incoming cosmic radiation. Below — 20 GeV/a a simulation was used by Ichimura et 
al. vi to compute the effect of the geomagnetic cutoff. Beyond ~ 20 GeV/a the predominant 
effect of decreasing energy resolution is to require a larger correction. There is still some 
sensitivity to the power index as shown v This is undesirable and is best 

controlled by making good measurements so that the energy resolution is small. The 
variation with the assumed power index of the primary radiation is relatively small if the 
value of log, 0 (E es /E, mc ) is in the range of 0.2 to 0.3 as was obtained by Ichimura et al." in 
their accelerator measurements. For similar cosmic ray species the relative flux ratios will 
be well determined but should the energy resolution change significantly with species 
both the relative and absolute intensities may not be reliable. 
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5. Conclusions 

The reduced opening angle technique offers a simple way to measure the charge 
and energy of cosmic rays. The technique has minimal model dependence at least up to 
several hundred GeV/a. As all the measurements would be made with the same optical 
techniques there are no energy normalization uncertainties. 

The energy range that can addressed is limited by the instrument dimensions. 1 he 
instrument can in principle be scaled to dimensions limited only by the lifting vehicle. In 
practice the needs for accurate geometric reference points is also a demanding constraint. 
For the geometry suggested here (flight path of 30 cm) energy measurements lo - 5 
TeV/a can be made with a charge resolution of —0.3 for the CNO group to ~ 1 for the iron 
group. Energy resolutions of -50% to + 80% have been obtained in accelerator 
calibrations at 200 GeV/a. A practical limit to the utility of the technique is set by the 
cosmic ray flux. The event rate per m 2 of the detector for iron cosmic rays exceeding 10 
TeV/a is 0.1 per day. 

The absolute flux measurements have a sensitivity to the assumed power index. 
The sensitivity is minimized if good energy resolution (log 10 (E £St /E tn J -0.2) is obtained. 

6. Appendices 


6.1 Appendix A: Electromagnetic dissociation 


Electromagnetic dissociation (EMD) occurs when the Lorentz contracted electric 
Held of the target nucleus contributes to the fragmentation of the projectile nucleus. In a 
simple picture EMD occurs when a virtual photon is exchanged between the target 
nucleus and the projectile often resulting in the excitation of a giant mullipole 
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resonance" 1 . The total cross-sections from electromagnetic processes can exceed die 
nuclear component but the charge changing electromagnetic cross-sections at least lor 
gold projectiles are not large as measured by Geer et al. xv ". The EMD component in such 
processes can be determined by the method of factorization where it has been shown that 
the partial a AZ (T,F) for a given projectile can be separated into factors y h depending only 
on the fragment and yT depending only on the target. 

<ru(T,F) = r y r T + <t imi ,(7\F) 


lb. 


where ct i;md (T,F) represents the EMD contribution. For this analysis the hydrogen partial 
cross-sections are used to define y r with yT=H=l . Hydrogen a are used since they should 
have minimal EMD. The other target a are divided by the hydrogen a to give. 




&IMI> 


(T, F) 


17 . 


To determine the magnitude ol 

yT a simple expression based on a fit to the cross-section is used. 


yT = 



where a„ (=2524 mb) is a normalization factor. Hence the EMD component can he 
determined by subtracting equation 17 from the total cross-section. Although m> 1 Ml) 
cross-sections for iron are reported by Geer et al. xv " we may estimate their magnitude as 

26 2 

lb. 
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Target 

Measured AZ = -1 
ct (mb) for l97 Au 
projectiles on a 
range of targets 

Measured AZ = -2 
a (mb) for 197 Au 
projectiles on a 
range of targets 


Deduced AZ = -1 
cr (mb) for 56 Fe 
projectiles on a 
range of targets 

Deduced AZ = -2 
ct (mb) for M 'Fe 
projectiles on a 
range of targets 

C 

15±10 

5±5 


1.65 

.55 

Al 

27±9 

11±5 


3.07 

1.21 

Cu 

55±10 

13±6 


6.05 

1.43 

Sn 

1 1 6± 1 2 

23±7 


12.76 

2.53 

Pb 

268±15 

42+7 


29.48 

4.62 


Table 7 Measured electromagnetic dissociation cross-sections for 197 Au projectiles 
and deduced values for iron (see text). 


From comparing Table 3 and Table 7 it is clear that the EMD contribution to the 
total charge changing cross-section is ~ 10 % for AZ = 1 and less for larger charge 
changes. For the opening angle technique the reaction channels with only one proton out 
produce essentially no deviation in the cosmic ray trajectory and no energy determination 
can be made. Hence we can neglect the EMD component to the reaction channels for this 
technique. 


6.2 Appendix B: Nuclear cross-sections and their energy dependence 
(E<10.6 GeV/a). 

In this appendix we discuss the energy dependence of the measured nuclear cross- 
sections and present some simple phenomenological models that have been fit to data. 
Nilsen et al. x have used a logarithmic scaling that was inspired by the logarithmic scaling 
used by the Particle Data Group™ 11 at high energies to describe nucleon/nucleon 
interactions. They parameterize the cross-sections as: 


<t(R p ,R t ) = xF(e\r p +R t -G(E) A/ff 


20 . 


F(E) = 1 + F t Ln(E ) + F 2 [Ln(E)] 2 

21 . 


G(E) = 1 + G, Ln(E) + G 2 [Ln(E)f 

22 . 


The F value scales the entire expression whilst the G value modifies the overlap 
term AR. Both terms depend on energy and on the model used for the overlap term AR. 
Two variants for AR are considered by Nilsen et al: equations 23 and 24. Fits are 
provided for each term and a series of measured conditions as shown in Table 8. These 
fits include data from Nilsen et al. (a) and from Binns et al. x,x (b) and Geer et al. xx (c ) 


A R=r 0 S 


23. 
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24. 


A R = 


rp5_ 

{a?+a?) 


Equation 

F, 

f 2 

G, 

g 2 

N 2 v 

25' 2 

-1.512±0.012 

1.11 ±0.04 

-1.4110.04 

0.9610.06 

4.65 

25 i3 

-1.61±0.16 

1.18+0.14 

-1.5110.19 

1.0410.18 

3.92 

25 14 

-0.28±0.10 

0.09±0.03 

-0.4010.10 

0.1110.04 

4.88 

26 12 

-1.24±0.19 

0.95±0.16 

-1.610.3 

1.110.2 

4.73 

26 13 

-1.28+0.18 

0.98±0.14 

-1.610.3 

1.110.2 

4.07 

26 14 ; 

-0.23±0.07 

0.09±0.02 

-0.6010.11 

0.2110.04 

4.86 


Table 8 Fit parameters for hard sphere models with overlap terms given by 
equations 23 and 24. 

The computed cross-sections have a strong energy dependence as shown in Table 
9. The necessity to use only the parameter specially developed for the particular energy 
regime of interest is clear. 


Equation 

Parameters 

Energy 

a barns 


None 

All 

4.509 

22 

25 13 

1.5 GeV/a 

4.69 

22 

25 13 

10.6 GeV/a 

6.115 

22 

26 14 

10.6 GeV/a 

2.701 


Table 9 Calculated cross-sections for iron projectiles on lead targets. Models with 
and with out energy dependence have been used - see text for details. 


6.3 Appendix C: Some measured fragment distributions 

In section 4.1 we present simple phenomenological formula for calculating 
nuclear fragment cross-sections. For the purposes of this study what may be of more 
interest are actual fragment yield data as presented by Waddington et al. XXI . Although 
these have been measured at low energies and are for ultra heavy nuclei the distributions 
are of interest. However, any conclusions drawn from them are subject to these caveats. 

In Figure 8 we show that the fragment yield (Z > 3) has little energy variation 
between 1 and 10.6 GeV/a. 


12 There are 43 values of Kr and Ag c(P,T,K) used in these fits from Nilsen et al’s. paper. 

13 These fits include the ten additional values of a(P,T,K) from Binns et al. 

14 These fits include the six additonal values at 10.6 GeV/a from Geer et al. with the electromagnetic 
contribution removed. 
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Figure 8 Comparison of the yield of Z > 3 
fragments at 1 and 10.6 GeV/a with a ,97 Au 
projectile. 



Figure 9 The fraction of a-particles to the 
total number of events at 1 and 10.6 GeV/a 
with a m Au projectile. 


In contrast to Figure 8 the yield of a-particles show some energy dependence. As 
the energy increases fewer events have no a-particles and there is a distinct increase in 
the region of 3-5 a-particles. The number of fragments plus a’s (Figure 10) shows a 
strong move towards less fragments as the energy increases and the yield of protons also 
increases with energy (Figure 11). 
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Figure 10 The yield of all fragments with 
Z > 2 as a function of fragment number 
for a m Au projectile. 



Figure 11 The ratio of the protons yield 
to the number of events as a function of 
the number of released protons for a 
197 Au projectile. 



pseudo rapidity’s of ~2.5 while the heavier 
minimum ionizing particles overlaps the Z > ^ 


It is also interesting to look at the 
pseudo rapidity (-logJtanG/2]) 

distributions which are shown in Figure 

12 . 

Figure 12 Pseudo rapidity 
distributions for a 197 Au projectile. 

Although the minimum ionizing 

particles (protons and pions) peak at 
fragments are at ~ 5.5 the tail of the 
! particles. Hence it is essential that the 
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emulsion films employed be able to unambiguously identify protons and pions from all 
other species. This is also a requirement to guard against wounded protons contaminating 
the distribution - see section 2.2 
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